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Abstract 

In contrast to hole-doped systems which have hole pockets centered at 
(±;r, in lightly electron-doped antiferromagnets the charged quasi- 

particles reside in momentum space pockets centered at (~,0) or (0,-). 
This has important consequences for the corresponding low-energy effec- 
tive field theory of magnons and electrons which is constructed in this 
paper. In particular, in contrast to the hole-doped case, the magnon- 
mediated forces between two electrons depend on the total momentum P 
of the pair. For P = the one-magnon exchange potential between two 
electrons at distance r is proportional to 1/r 4 , while in the hole case it 
has a 1/r 2 dependence. The effective theory predicts that spiral phases 
are absent in electron-doped antiferromagnets. 
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1 Introduction 



Although they are not yet high-temperature superconductors, understanding lightly 
doped antiferromagnets is a great challenge in condensed matter physics. A lot is 
known about hole- and electron-doped systems both from experiments and from stud 
ies of microscopic Hubbard or i-J-type models El B, I E H, M B , B , [j[ [j| 




Based on the work of Haldane [2J and of Chakravarty, Halperin, and Nelson [6j who 
described the low-energy magnon physics by a (2 + l)-d 0(3)-invariant nonlinear o- 
model, several attempts have been made to include charge carriers in the effective 
theory 0, 0, [hJ 12 1 . However, conflicting results have been obtained. For example, 
the various approaches differ in the fermion field content of the effective theory and 
in how various symmetries are realized on those fields. In particular, it has not yet 
been established that any of the effective theories proposed so far indeed correctly de- 
scribes the low-energy physics of the underlying microscopic systems in a quantitative 



manner. 



In analo gy t o chiral perturbation theory for the pseudo-Nambu-Goldstone pions 
of QCD [34], |35j, the (2 + l)-d 0(3)-invariant nonlinear cr-model has been established 
as a systematic and quantitati vely correct low-energy effective field theory in the 



pure magnon sector 
perturbation theory [4 




37 
46 
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47 




In analogy to baryon chiral 
the effective theory for pions and nucleons 



we have recently extended the pure magnon effective theory by including charge 



carriers 



49, 50, 51 



The effective theory provides a powerful theoretical framework 
in which the low-energy physics of magnons and charge carriers can be addressed in a 
systematic manner. The predictions of the effective theory are universal and apply to a 
large class of doped antiferromagnets. This is in contrast to calculations in microscopic 
models which usually suffer from uncontrolled approximations and are limited to just 
one underlying system. While some results obtained with the effective theory can be 
obtained directly from microscopic systems, the effective field theory treatment allows 
us to derive such results in a systematic and more transparent manner and it puts them 
on a solid theoretical basis. In order not to obscure the basic physics of magnons and 
charge carriers, the effective theory has been based on microscopic systems that share 
the symmetries of Hubbard or t-J-type models. In particular, effects of impurities, 
long-range Coulomb forces, anisotropies, or small couplings between different Cu02 
layers have so far been neglected, but can be added whenever this becomes desirable. 
Before such effects have been included, one should be aware of the fact that the 
effective theory does not describe the actual materials in all details. Still, for systems 
that share the symmetries of the Hubbard or t-J model, the effective theory makes 
predictions that are exact, order by order in a systematic low-energy expansion. 

Hole-doped cuprates have hole pockets centered at lattice momenta (±^,±^). 
The location of the hole pockets has important consequences for the fermion field 
content of the effective theory and on the realization of the various symmetries of these 
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Figure 1: The dispersion relation E(p) of a single electron in the t-t'-J model (on a 
32 x 32 lattice for J = OAt and t' = —0.3t) with electron pockets centered at (-,0) 
and (0,f). 



fields. In electron-doped cuprates the charged q uasi par ticles reside in momentum 
space pockets centered at (f, 0) or (0, f ) 0, 0, 0, H l3ll . |32j, |33| . We have computed 
the single-electron dispersion relation in the t-t'-J model shown in figure 1. The energy 
E(p) of an electron is indeed minimal when its lattice momentum p = (pi,p2) is located 
in an electron pocket centered at (-, 0) or (0, -). The location of these pockets again 
has important effects on the electron dynamics, which turns out to be quite different 
from that of the holes. In particular, in contrast to hole-doped systems, in electron- 
doped antiferromagnets the magnon-mediated forces between two electrons depend 
on the total momentum P of the pair. For P = the one-magnon exchange potential 
between two electrons at distance r is proportional to 1/r 4 , while in the hole case it 
has a 1/r 2 dependence. The different locations of electron and hole pockets also affect 
the phase structure. While spiral phases are possible in the hole-doped case 

BBS, 



they are absent in electron-doped cuprates [17], |20j, [22j, [24j, [26 



The paper is organized as follows. In section 2 the symmetries of charge carrier 
fields are summarized. Based on this, the electron fields are identified and the hole 
fields are eliminated. The low-energy effective action for magnons and electrons is 
then constructed in a systematic manner. Section 3 contains the derivation of the 
one-magnon exchange potential between two electrons as well as a discussion of the 
corresponding Schrodinger equation. In section 4 spiral configurations of the staggered 
magnetization and in section 5 the reduction of the staggered magnetization upon 
doping are investigated. Section 6 contains our conclusions. The somewhat subtle 
transformation of the one-magnon exchange potential from momentum to coordinate 
space is discussed in an appendix. 
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2 Symmetries of Magnon and Electron Fields 



In this section, based on [49j, [51[, we summarize the transformation properties of 
magnon and charge carrier fields. We then identify the electron fields and eliminate 
the hole fields in order to construct the low-energy effective theory for magnons and 
electrons. 



2.1 Symmetries of Magnon Fields 



In an antiferromagnet the spontaneous breaking of the SU(2) S spin symmetry down to 
U(l) s gives rise to two massless magnons. The staggered magnetization is described 
by a unit- vector field 

e(x) = (e\(x) , e2(x) , es(x)) = (sin#(x) cos ip(x), sin 9(x) sin ip(x), cos8(x)), (2.1) 

in the coset space SU(2) S /U(1) S = S 2 , where x = (xi,x 2 ,t) is a point in (2 + l)-d 
space-time. It is convenient to use a CP(1) representation in terms of 2 x 2 Hermitean 
projection matrices P(x) with 



P{x) = -(1 + e(x) ■ a), P(x) j = P(x), TiP(x) = 1, P(x) 2 = P(x) 



(2.2) 



As discussed in detail in [49[, the magnon field transforms as 
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(2.3) 



The various symmetries are the SU(2) S spin rotations, the non-Abelian SU(2)q exten- 
sion of the U(1)q fermion number symmetry (also known as pseudo-spin symmetry) 
that arises in the Hubbard model at half-filling, the displacement symmetry by one 
lattice spacing in the i-direction Di, the symmetry Di combined with the spin rotation 
io 2 resulting in D\, as well as the 90 degrees rotation O, the reflection at the Xi-axis 
R, time reversal T, and T combined with the spin rotation ia 2 resulting in T' . 

The spontaneously broken SU(2) S symmetry is nonlinearly realized on the charge 
carrier fields. The global SU(2) S symmetry then manifests itself as a local U(l) s sym- 
metry in the unbroken subgroup, and the charge carrier fields couple to the magnon 
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field via composite vector fields. In order to construct these vector fields one first 
diagonalizes P(x) by a unitary transformation u(x) G SU(2), i.e. 



u(x)P(x)u(xY = -(1 + (T 3 ) = 



1 




tin (a;) > 0, 



u(x) 



cos ^2 sin m eX p(-^(x)) 



cos 



e(x) 



(2.4) 



2 

— sin ^ exp(i(p(x)) ,, 
Under a global SU(2) S transformation g, the diagonalizing field u(x) transforms as 

u(x)' = h(x)u(x)g\ Unix)' > 0. (2.5) 
This defines the nonlinear symmetry transformation 



h(x) = exp(ia(x)cr 3 ) 



exp(ia(x)) 






exp(— ia(x)) 



e U{i) t 



(2.6) 



Under the displacement symmetry Di the staggered magnetization changes sign, i.e. 
Di e(x) = —e(x), and one obtains 



u(x) = t(x)u(x), t(x) = 



— exp(— i(p(x)) 

exp(i(p(x)) 



(2.7) 



In order to couple magnons and charge carriers, one constructs the traceless anti- 
Hermitean field 

v^(x) = u(x)d ll u(x)\ (2.8) 

which transforms as 



SU(2) S 




= h{x)[v^{x) + d fl ]h(x)\ 




SU(2) Q 


®v^x) 


= Vn(x), 




Di 


Dl v^(x) 


= t(x)[v^(x) + <^]r(x) t , 




D[ 


D' 1 \ 


= v^x)*, 







°Vi(x) 


= EijVjtOx), °v t (x) = 


v t (Ox), 


R 


R v 1 (x) 


= vi{Rx), R v 2 {x) = 


-v 2 (Rx), R v t (x) 


T 


T Vj (x) 


= Di Vj (Tx), T v t (x) = 


- Di v t (Tx), 


T' 


T' 1 \ 
Vj{x) 


= D i Vj (Tx), T 'v t (x) = 


- D *v t (Tx). 



(2.9) 

The field v^x) decomposes into an Abelian "gauge" field v^(x) and two "charged" 
vector fields v^(x), i.e. 



v^x) = ivfa)a a , vfa) = vfa) =f ivf^x). 



(2.10) 
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Figure 2: Eight lattice momenta and their periodic copies. In the cuprates the holes 
reside in momentum space pockets centered at lattice momenta (if-, if 1 ) which are 
represented by the four crosses, while electrons reside at ,0) or (0, ^) (represented 
by the circles). 

2.2 Fermion Fields in Momentum Space Pockets 



In [5l| matrix-valued charge carrier fields 



i/j_(x) 



4)Z h '(x) -?p+ k \x) 



(2.11) 



have been constructed. Here k! = k + (-, -) and i/j±(x) and ip±(x) are independent 
Grassmann fields which are associated with the following eight lattice momentum 
values illustrated in figure 2 



*=(*.,*,)€ {(0,0), (-,-), (-„), (0,-),( ± -, ± -)}. 
The charge carrier fields transform as 



;2.i2) 



SU(2) S : 


tf fc (z)' 


SU(2) Q : 




A: 




^: 




: 




i2 : 


R ^ k (x) 


T : 


T ^\x) 




T ^ k Hx) 



h{x)V k {x), 
^ k (x)tt T , 

exp(ikia)r(x)^ k (x)a3, 
eiq)(ik i a)(ia 2 )'^ k (x)a 3 , 
^ 0k (Ox) } 
V Rk {Rx), 

r{Tx){ia 2 ) [^- k \Tx) T ] a 3 , 
-a 3 [V- k (Tx) T ] (ia 2 )V(Tx 
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T : T V(x) = - [m- k \Tx) T ] a 3 , 

T 'm k \x) = a 3 [^ k (Tx) T ] . (2.13) 

Here Q G SU(2)q and Ok and Rk are the momenta obtained by rotating or reflecting 
the momentum k. 



2.3 Electron Field Identification and Hole Field Elimination 



ARPES measurements as well as theoretical investigations [271, l2_8|, |29|, [3CJ, [3JJ, |32|, [33 
(see also figure 1) indicate that electrons doped into an antiferromagnet appear in 
momentum space pockets centered at 



fc=(-,0), k'=(0,-). (2.14) 
a a 



Hence, only the fermion fields with these two momentum labels will appear in the low- 
energy effective theory. Using the transformation rules of eq. fl2.13p one can construct 
the following invariant mass terms 

~Tr[A4(* fet <7 3 * fe ' + ^ fc 'V 3 ^ fc ) + m{^ k ^ k a 3 + $ k '^ k 'a 3 

+ ^)(^). (2,5) 

The terms proportional to are S , f/(2)Q-invariant while those proportional to m are 
only [/(l)g-invariant. By diagonalizing the mass matrices, electron and hole fields 
can be identified. The resulting eigenvalues are m ± M.. In the Si7(2)Q-symmetric 
case, i.e. for m = 0, there is an electron-hole symmetry. The electrons correspond 
to positive energy states with eigenvalue M. and the holes correspond to negative 
energy states with eigenvalue — M.. In the presence of S77(2)Q-breaking terms these 
energies are shifted and electrons now correspond to states with eigenvalue m + A4, 
while holes correspond to states with eigenvalue m — Ai. The electron fields are given 
by the corresponding eigenvectors 

^+(ar) = ^ +1>+(x)], t-W = ^[^(x)-^'W]. (2.16) 



7 



Under the various symmetries they transform as 



OT 7Yo\ 

bU(2) s 


ip±{x) 


= exp{±ta{x))ip±{x), 


U(1)q 


i I \ 

w ip±{x) 


= exp(ujL))ip±{x), 


A 




= =F exp(i/cja) exp(^fiLp(x))ip T (x), 


n' 

D i 




= ± exp(ikia)^pzf(x), 


r\ 
yj 


ip±{x) 


l_„/. /'n^'i 


t? 
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lp±{X) 


— ip±{UX), 


1 


i>±{x) 


= exp(^np(Tx))ip ] ± {Tx), 






= -exp(±i(p(Tx))if>±(Tx), 


T' 




= -4>l(Tx), 






= ^±{Tx). 



The action of magnons and electrons must be invariant under these symmetries. 



2.4 Effective Action for Magnons and Electrons 

We decompose the action into terms containing different numbers of fermion fields 
(with even) such that 

S[rpL^±,P] = J d 2 x dt (2.18) 

The leading terms in the effective Lagrangian without fermion fields are given by 

C = p s Tr [d t Pd,P + ^d t Pd t P] , (2. 19) 

with the spin stiffness p s and the spinwave velocity c as low-energy parameters. The 
terms with two fermion fields (containing at most one temporal or two spatial deriva- 
tives) describe the propagation of electrons as well as their couplings to magnons, and 
are given by 



s=+, 



+ iK{D 1 ^ s v s l ^. s - il)\v{D^_ s - D^\v s ^. s + ij\v s 2 D 2 ^ s )j . (2.20) 

Here M is the rest mass and M' is the kinetic mass of an electron, K is an electron- 
one-magnon, and N is an electron-two-magnon coupling, which all take real values. 
The covariant derivatives are given by 

D t ip±(x) = [d t ± ivt(x) - fi\ ip±(x), 

Diip±(x) = [di ± ivf(x)] rl>±(x). (2.21) 
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The chemical potential /x enters the covariant time- derivative like an imaginary con- 
stant vector potential for the fermion number symmetry U(1)q. 



Next we list the contributions with four fermion fields including up to one temporal 
or two spatial derivatives 

s=+- „ 

G 

+ ^(Drftyrflrf-, - tiDtMtrf-s)] ■ (2.22) 

Since it contains D t , the term proportional to G$ would imply a deviation from canon- 
ical anticommutation relations in a Hamiltonian formulation of the theory. Fortu- 
nately, this term can be eliminated by a field redefinition ip s — > ip s + ^■ip s ip_ s il>^ s . 
The redefined field obeys the same symmetry transformations as the original one and 
is constructed such that after the field redefinition Gs = 0. All other terms in the 
action are reproduced in their present form. 

For completeness, we finally list the only contribution with more than four fermion 
fields, again including up to one temporal or two spatial derivatives 

£ 6 = HDtflDiMlMtrf.,. (2.23) 
«=+,- 

The leading fermion contact term is proportional to G±. Due to the large number of 
low-energy parameters G 2 , ■■■,Gf, H, the higher-order terms are unlikely to be used 
in practical applications. We have used the algebraic program FORM [53j], and inde- 
pendently thereof, the GiNaC framework for symbolic computation within the C++ 
programming language 54], to verify that the terms listed above form a complete 
linearly independent set. 

It should be noted that, unlike in the hole case, the leading terms in the effective 
action are not invariant against Galilean boosts. This is not unexpected because the 
underlying microscopic systems also lack this symmetry. The lack of Galilean boost 
invariance has important physical consequences. In particular, the magnon-mediated 
forces between two electrons will turn out to depend on the total momentum P of 
the pair. Thus it is not sufficient to consider the two particles in their rest frame, 
i.e. at P = 0. This is due to the underlying crystal lattice which defines a preferred 
rest-frame (a condensed matter "ether"). 
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3 Magnon- mediated Binding between Electrons 



We treat the forces between two electrons in the same way as the ones in the effec- 



tive theory for magnons and holes |50l . l5l| . As in that case, one- magnon exchange 



dominates the long-range forces. In this section we calculate the one-magnon ex- 
change potential between two electrons and we solve the corresponding two-particle 
Schrodinger equation. 



3.1 One-Magnon Exchange Potential between Electrons 

In order to calculate the one-magnon exchange potential between two electrons, we 
expand in the magnon fluctuations m^x), m 2 (x) around the ordered staggered mag- 
netization, i.e. 



e>) = (^,^,lWo(m 2 ) 



■J^(x) = — dp [m 2 (x) ± influx)] + C(m 3 ), 

Ps 



v^(x) = — [m 1 (x)d^m 2 (x) - m 2 (x)d fl m 1 (x)] + C(m 4 ). (3.1) 

The vertices with v^(x) (contained in D^) involve at least two magnon fields. Hence, 
one-magnon exchange results exclusively from vertices with v^(x). Thus, two electrons 
can exchange a single magnon only if they have antiparallel spins (+ and — ), which 
are both flipped in the magnon exchange process. We denote the momenta of the 
incoming and outgoing electrons by p± and p± , respectively. Furthermore, q represents 
the momentum of the exchanged magnon. We also introduce the total momentum P 
as well as the incoming and outgoing relative momenta p and p 1 

P = p + + P- = p+ + P-, 

P=^(P+~P~), P'=l(p+-P-)- (3-2) 

Due to momentum conservation we then have 

q = p + p'. (3.3) 

Figure 3 shows the Feynman diagram describing one-magnon exchange. In momentum 
space the resulting one-magnon exchange potential takes the form 

K 2 1 

(P+ P-\V\P+P-) = 7> 2 fa? ~ ?2 + 2 (<?iP-i - 92P-2)] [?? -<&- 2(giP+i - g 2 P+2)] 

^Ps Q 

xd(p + +p^-p;-pl). (3.4) 
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Figure 3: Feynman diagram for one-magnon exchange between two electrons with 
antiparallel spins undergoing a spin-flip. 

Transforming the potential to coordinate space is not entirely trivial and is thus 
discussed in the appendix. In coordinate space the resulting potential is given by 



(r\rL\V\r + r_) 



K 



2 



2-Kp s 



l2 cos(4^) P 2 cos(2(y + x)) 
r 4 2 r 2 



5{r + -rL) 8{r„-rl). (3.5) 



Here p is the angle between the distance vector f = f + — r_ of the two electrons and 
the Xi-axis. In contrast to the hole case, the potential depends on the magnitude 
P of the total momentum P as well as on the angle \ between P and the xj-axis. 
For P = the one-magnon exchange potential between two electrons falls off as 
1/r 4 , while in the hole case it is proportional to 1/r 2 . Retardation effects enter at 
higher orders only and thus the potential is instantaneous. We have omitted short- 
distance 5-function contributions to the potential which add to the 4-fermion contact 
interactions. Since we will model the short- distance repulsion by a hard-core radius, 
the ^-function contributions will not be needed in the following. 



3.2 Schrodinger Equation for two Electrons 



Let us consider two electrons with opposite spins + and — . The wave function depends 
on the relative distance vector r which points from the spin — electron to the spin 
+ electron. Magnon exchange is accompanied by a spin-flip. Hence, the vector r 
changes its direction in the magnon exchange process. The resulting Schrodinger 
equation then takes the form 



1 K 2 
A*(r)H 
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cos(4y?) P 2 cos(2(y? + x)) 





\e p2 l 


\p(-f) = 






2M'_ 



*(r). (3.6) 



For simplicity, instead of explicitly using the 4-fermion contact interactions, we model 
the short- distance repulsion between the electrons by a hard-core of radius r , i.e. we 
require \l/(f) = for \r\ < To- In contrast to the hole case 0, 51], we have not been 
able to solve the above Schrodinger equation analytically. Instead, we have solved it 
numerically. A typical probability distribution for the ground state is illustrated in 
figure 4 for P = 0. The probability distribution resembles d xy symmetry. However, 
due to the 90 degrees rotation symmetry the continuum classification scheme of an- 
gular momenta is inappropriate. Under the group of discrete rotations and reflections 
the ground state wave function transforms in the trivial representation. 
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Figure 4: Probability distribution for the ground state of two electrons with total mo- 
mentum P = (0, 0). 

Due to the lack of Galilean boost invariance, the two-electron bound state changes 
its structure when it is boosted out of its rest frame. Of course, an electron pair with 
total momentum P^O costs additional kinetic energy P 2 /2M' for the center of mass 

— * 

motion. In addition, the binding energy also depends on P. The strongest binding 

— * 

arises when the total momentum P points along a lattice diagonal. The corresponding 
probability distribution is illustrated in figure 5. Since they depend crucially on the 
precise values of the low-energy parameters, we have not attempted an extensive 
numerical investigation of the binding energy and other properties of the two-electron 
bound states. Once the low-energy parameters have been determined for a concrete 
underlying microscopic system, a precise calculation of the physical properties of the 
two-electron bound state is straightforward using the numerical method employed 
above. 

In order to gain at least some approximate analytic insight into the bound state 
problem, let us also consider the semi-classical Bohr-Sommerfeld quantization. First, 
we consider a pair of electrons with total momentum P = moving relative to each 
other along a lattice diagonal. The classical energy of the periodic relative motion is 
given by 

E = M>(Pl-**<. (3.7) 
\dt) np s r A 1 ; 
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Figure 5: Probability distribution for the ground state of two electrons with total mo- 
mentum P = 4=(P, P) along a lattice diagonal. 



The Bohr-Sommerfeld quantization condition implies 

2 n ts2 



S + ET 







[ dt 




Jo 





dr 
dt 



6K 2 
np s r 4 



E 



T 

2 I dt W 
'o 



dr 

~dl 



if drM'%=\\ dr J. 

Jro dt Jro V 



EM 1 + 



6K 2 M' 
np s r A 



27rn, 



(3.8) 



where S is the action, T is the period of the motion, and n is a positive integer. The 
hard-core radius r is a classical turning point and R is the other classical turning 
point determined by 

* = (3-9) 

The above equations lead to a relatively complicated expression for the energy in terms 
of elliptic integrals. Instead of investigating these expressions, we limit ourselves to 
estimating the number of bound states. For this purpose, we set E = which implies 
R = oo, and we then obtain 



n 



dr 



ro 



I24:K 2 M' 
Ti 3 p s r 4 



! 24K 2 M' 
n 3 p s r 2 



(3.10) 



The brackets denote the nearest integer smaller than the expression enclosed in the 
brackets. In particular, Bohr-Sommerfeld quantization suggests that a bound state 
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exists only if 



2AK 2 M' 
ir 3 p s r 2 



> 1. 



(3.11) 



Of course, one should be aware of the fact that this is at best a semi-quantitative 
estimate because Bohr-Sommerfeld quantization should not be trusted quantitatively 
for small quantum numbers. Let us also repeat these considerations for Again, 
we consider P = 4j(P, P) such that the diagonal motion of an electron-pair has the 
energy 

'dr\ 2 K 2 (\2 P 2 
7 + 



E — M' 



dt 



2np s 



2r 2 



(3.12) 



In complete analogy to the P = case one then obtains 



n 



dr 



ro 



'2K 2 M' 
n 3 p s 



12 P 2 
~^ + 2r 2 



OO, 



(3.13) 



which suggests that infinitely many two-electron bound states exist for P ^ 0. This 
is similar to the two- hole problem which has a 1/r 2 potential with infinitely many 
bound states already for P = ^0, 51 . 
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have been considered before by Kuchiev 
in the context of the t-t'-J model. In contrast to the hole-case 13 



Two-electron bound states with P 
and Sushkov 



with a 1/r 2 potential and an infinite number of bound states, in the electron-case 
only a finite number of bound states was found. While some results of that work 
agree qualitatively with the results of our effective theory, there are also important 
differences. For example, in 19j the magnon-electron vertex was considered to be the 
same as the magnon-hole vertex, while the two vertices are different in the effective 
theory. 



4 Investigation of Spiral Phases 

In the following we will investigate phases with constant fermion density. The most 
general magnon field configuration e(x) which provides a constant background field 
for the doped electrons is not necessarily constant itself, but may represent a spiral in 
the staggered magnetization. While a spiral costs magnetic energy proportional to the 
spin stiffness p s , the electrons might lower their energy by propagating in the spiral 
background. However, we will find that spiral phases are not energetically favorable 
in electron-doped systems. 
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4.1 Spirals with Uniform Composite Vector Fields 



Since the electrons couple to the composite vector field Vi(x) in a gauge covariant 
way, in order to provide a constant background field for the electrons, Vi(x) must be 
constant up to a gauge transformation, i.e. 



vf(x)' = vf(x) — dia(x) = diip(x) sin 



. 2 e(x) 



dja(x) = cf. 



vf(x)' = vf(x) exp ( ± 2ia(x)) 



- [diip(x) sin9(x) ± idi6(x)\ exp ( ± 2ia(x) =F ipfe)) = c t 



(4.1) 



with cf and cf being constant. As shown in 52|, the most general configuration 
that leads to a constant Vi(x)' represents a spiral in the staggered magnetization. In 
addition, by an appropriate gauge transformation one can always put 



cf = c- 



The magnetic energy density of such configurations takes the form 
e m = ^-die(x) ■ die(x) = 2p s vf(x)v~(x) = 2p s dCi. 



(4.2) 
(4.3) 



We now consider a concrete family of spiral configurations with 

9{x) = O , f(x) = kiXi, 

which implies 



Q 

v t (x) = 0, vf(x) = ki sin 2 — 
Performing the gauge transformation 



k 

v ii x ) = 77 si 11 ex P(T^i^i)- 



(4.4) 



(4.5) 



one arrives at 



such that 



v t {x)' = v t (x) - d t a(x) = 0, 
vf(x) — dja(x) = h 



vfixY 



2 9 1 



vf(x)' 



sm 

2 2 

k; 



vf{x) exp ( ± 2ia(x)j = sin^o = c. 



cos 6n 



— tan^n- 



The magnetic energy density then takes the form 

e m = 2p s dCi = ^(k\ + kf) sin 2 O . 



(4.6) 



(4.7) 
(4.8) 

(4.9) 
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4.2 Fermionic Contributions to the Energy 



Let us now compute the fermionic contribution to the energy, first keeping the param- 
eters cf and Q of the spiral fixed, and neglecting the 4-fermion contact interactions. 
The Euclidean action of eq. (12.201) implies the following fermion Hamiltonian 



L 2M' 

s=+,- 

+ ^(D^tf _ s _ *t v j Dl ^_ a _ D 2 *t v ^_ 8 + ^^^j ? (4 . 10) 

with the covariant derivative 

Di* ± (x) = [di ± ^ 3 (x)]^±(x). (4.11) 

Here *&±{x) and \l/±(x) are creation and annihilation operators (not Grassmann num- 
bers) for electrons with spin parallel (+) or antiparallel (— ) to the local staggered 
magnetization. The Hamiltonian is invariant under time-independent U(l) s gauge 
transformations 

W±{x)' = exp ( ± ia(x))^ ± (x), 
vf(x)' = vf(x) — dia(x), 



exp(±2ia(x)). (4.12) 

We now consider electrons propagating in the background of a spiral in the staggered 
magnetization with 

vUx)' = cl vt(x)' = Cl eR. (4.13) 

After an appropriate gauge transformation, the fermions propagate in a constant 
composite vector field Vi{x)' . The Hamiltonian is diagonalized by going to momentum 
space. The Hamiltonian for single electrons with spatial momentum p = {p\,P2) is 
given by 

„3>,2 



M + + Nc i°i 2K (-PlCl + V2C2 

2K(-p 1 c 1 + p 2 c 2 ) 
The diagonalization of the Hamiltonian yields 



H(p)=\^' 2M' ' — \ 1 3^2 ' ^ wy . (4.14) 
2K(-p lCl +p 2 c 2 ) M+&£U-+Nc i c i 



2 / 3 \ 2 / / 3 \ ^ 

E±(p) = M + Pi ^ i} + Nc lC% ± J (^j + 4Ki( PlCl - p 2 c 2 f. (4.15) 

It should be noted that in this case the index ± does not refer to the spin orientation. 
In fact, the eigenvectors corresponding to E±(p) are linear combinations of both spins. 
Since cf does not affect the magnetic contribution to the energy density, it can be fixed 
by minimizing E±(p) which leads to cf = c\ = 0. According to eq. fl4.7p this implies 
that 8 = |, i.e. the spiral is along a great circle on the sphere S 2 . For cf = the 
energies of eq. fl4.15p reduce to 

2 

E ± (p) = M+^ + Nc lCi ± 2K\ PlCl - p 2 c 2 \. (4.16) 
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Figure 6: Lines of constant energy for electrons propagating in a spiral configuration. 
The contours of the lower energy E_ (p) are shown on the left, and the contours of the 
higher energy E + {p) are displayed on the right. 



The lines of constant energy are shown in figure 6. In particular, the lines of constant 
E_{p) are circles centered around ±2KM'(ci, — c 2 ). For given q ^ we now fill the 
lowest energy states with a small number of electrons. The filled electron pockets 
are circles centered around ±2KM'(ci, — c 2 ) with a radius determined by the kinetic 
energy 

T = 2W< T 2KM ' C ^ + ^ 2 ± 2KM ' C 2) 2 ] (4-17) 
of an electron at the Fermi surface. The two occupied circular electron pockets define 
a region P in momentum space. The area of this region determines the fermion density 

as 



1 f ,0 1 



n 



The two circles do not overlap as long as n < ^M' 2 K 2 CiCi. The kinetic energy density 
of the filled region P is given by 

1 



t = 



(4.19) 



and the total energy density of electrons then is 



77-77 ^ 

e e = (M + Naci - 2K 2 M'c t c t )n + —. (4.20) 

The resulting total energy density that includes the vacuum energy density e as well 
as the magnetic energy density e m is given by 



im 



2 



e = e + e m + e e = e + 2p sCi c t + (M + Naa - 2K 2 M' Cl c t )n + — -. (4.21) 



2M 

17 




Figure 7: The energies E_(p) (solid curve) and E + (p) (dashed curve) along the line 
p oc (ci, — c 2 ) (the dashed lines in figure 6) define two independent parabolic dispersion 
relations. 



For p s > (K 2 M' — \N)n (which is always satisfied for sufficiently small density n) the 
energy is minimized for q = and the value of the energy density at the minimum is 
given by 

7T77 2 

e = e + Mn + — - . (4.22) 



However, one should not forget that, as the q become smaller, the two occupied 
circles eventually touch each other once -M' 2 K 2 CiCi = n. Interestingly, in this mo- 
ment the states with energy E + (p) also become occupied. Indeed, as one can see 
in figure 6, the almond-shaped region of occupied states with energy E + {p) and the 
peanut-shaped region of occupied states with energy E_(p), combine to two complete 
overlapping circles. This is also illustrated in figure 7, in which the energies E_(p) 
and E + (p) combine to form two overlapping parabolic dispersion relations. As a re- 
sult, eq. fl4.22p is still valid even when the occupied circles overlap. Consequently, the 
energy minimum is indeed at q = and thus a homogeneous phase arises. This is 
in contrast to hole-doped cuprates for which a spiral phase is energetically favored 



for intermediate values of the spin stiffness p s 52] . The effective theory predicts that 



spiral phases are absent in electron-doped antiferromagnets. 



4.3 Inclusion of 4-Fermion Couplings 

Let us also calculate the effect of the 4-fermion contact interactions on the energy 
density. We perform this calculation to first order of perturbation theory, assuming 
that the 4-fermion interactions are weak. Depending on the underlying microscopic 
system such as the Hubbard model, the 4-fermion couplings may or may not be small. 
We like to point out that, while the on-site Coulomb repulsion responsible for anti- 
ferromagnetism is always large in the microscopic systems, the 4-fermion couplings in 
the effective theory may still be small. If they are large, the result of the perturbative 
calculation should not be trusted. 
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The perturbation of the Hamiltonian due to the leading 4-fermion contact term of 
eq. fl2.22p is given by 

AH = I d 2 x ^s^-s. (4.23) 

s=+- 

It should be noted that ^l(x) and ^ s (x) again are fermion creation and annihilation 
operators (and not Grassmann numbers). The terms proportional to G2, G3, G-j 
are of higher order and will hence not be taken into account. The fermion density is 
equally distributed among the two spin orientations such that 

= = ~. (4.24) 

The brackets denote expectation values in the unperturbed state determined before. 
Since the fermions are uncorrelated we have 

(*l* a *l a *- 8 ) = <*L.tt_ a ). (4.25) 

Taking the 4-fermion contact terms into account in first order perturbation theory, 
the total energy density of eq. ( 14.221) receives an additional contribution and now reads 

e = 6o + Mn+ — + -V. (4.26) 



5 Reduction of the Staggered Magnetization upon 
Doping 

The order parameter of the undoped antiferromagnet is the local staggered magne- 
tization M a (x) = M. s e(x) with M. s being the length of the staggered magnetization 
vector. In a doped antiferromagnet the staggered magnetization receives additional 
contributions from the electrons such that 



M s (x) 



M s — m \J ^Kx)ip s (x) e(x). (5.1) 



The low-energy parameter m determines the reduction of the staggered magnetization 
upon doping. Further contributions to M s (x) which include derivatives or contain 
more than two fermion fields are of higher order and have thus been neglected. Using 

s=+- 

we then obtain 

M s (n) = M s - mn, (5.3) 

i.e. at leading order the staggered magnetization decreases linearly with increasing 
electron density. The higher-order terms that we have neglected will give rise to 
sub-leading corrections of 0(n 2 ). 
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6 Conclusions 



In analogy to the hole-doped case [49|, |5Jj, we have constructed a systematic effective 
field theory for lightly electron-doped antiferromagnets. Interestingly, the different 
locations of electron- and hole-pockets in the Brillouin zone have important conse- 
quences for the dynamics. 

In the hole-doped case, the pockets are located at (±^, which gives rise to 
a flavor index that determines to which pocket a hole belongs. Due to spontaneous 
symmetry breaking, holes and magnons are derivatively coupled. The leading magnon- 
hole coupling contains a single spatial derivative and is responsible for a variety of 
interesting effects. First, it leads to a 1/r 2 potential between a pair of holes which 
gives rise to an infinite number of two-hole bound states 5(3, 51]. Remarkably, in the 



hole-doped case, in the c — > oo limit the symmetries give rise to an accidental Galilean 
boost invariance. Hence, it is sufficient to consider the bound state in its rest-frame. 
Second, in the hole-doped systems, the single-derivative magnon-hole coupling gives 
rise to a spiral phase for intermediate values of p s . 

In the electron-doped case discussed in this paper, the momentum space pockets 
are located at (~,0) and (0, ^). Due to antiferromagnetism these are actually two 
half-pockets which combine to a single electron-pocket. Hence, in contrast to the 
hole-doped systems, electrons do not carry an additional flavor index. As in the hole- 
case, electrons and magnons are derivatively coupled. However, due to the different 
implementation of the symmetries, the leading magnon-electron coupling now con- 
tains two spatial derivatives. In other words, at low energies magnons are coupled to 
holes more strongly than to electrons. As a consequence, the one-magnon exchange 
potential between two electrons in their rest-frame decays as 1/r 4 and is hence weaker 
at large distances than in the hole-case. Still, magnon exchange is capable of binding 
electrons. As another consequence of symmetry considerations, an accidental Galilean 
boost invariance is absent in the electron-case. Indeed, the one-magnon exchange po- 
tential depends on the total momentum P of the electron-pair, and it is hence not 
sufficient to consider the system in its rest-frame. The momentum-dependent con- 
tribution to the potential is proportional to P 2 /r 2 , which gives rise to a non-trivial 
structure of moving bound states. As another consequence of the weakness of the 
magnon-electron coupling, in contrast to the hole-doped case, spiral phases are ener- 
getically unfavorable for electron-doped systems. While this is not a new result, we 
find it remarkable that it follows unambiguously from the very few basic assumptions 
of the systematic low-energy effective field theory, such as locality, symmetry, and 
unitarity. 

We like to point out that the systematic effective field theory approach is univer- 
sally applicable to a large class of antiferromagnets. While it remains to be seen if the 
effective theory can also be applied to high-temperature superconductors, it makes 
unbiased, quantitative predictions for both lightly hole- and electron-doped cuprates 
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and should be pursued further. 
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A Magnon Exchange Potential in Coordinate 
Space 

In this appendix we discuss the transformation of the one-magnon exchange potential 
between two electrons from momentum space to coordinate space, which is not entirely 
straightforward. 

In momentum space the one-magnon exchange potential is given by 



V{p,p') = [q{ - q\ + 2(q 1 p_ 1 - q 2 p- 2 )] [q\ -<&- - ^2^+2)] • (A.2) 



{p+p-\y\p+p-} = v(p,p')5(p + +p_ -pi-p-), 



(A.l) 



with 



Using 



P— * 1 — * — * / 1 — * / 
= P++P- = P+ + P-, 



P = o 




P' = 7,(P+-P-), 



q = p + p', 



(A.3) 



it is easy to show that 



V(p,p') = V (p,p') + Vp(p,p'), 



(A.4) 



with the rest-frame potential 



Vo(p,p') 



K 2 pj-Pi -pl+p' 2 2 

IPs (Pi + Pi)' + (P2 + P' 2 ) 



(A.5) 



and the momentum- dependent contribution 



Vp(p,p') 



K 2 [Pi(pi+p / 1 )-A(P2+P , 2 )] 
2p s (Pi+P'i) 2 + (P2 + P' 2 ) 2 



(A.6) 
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The potential in coordinate space is the Fourier transform of the potential in momen- 
tum space 

V(x,x') — -. — — [ d 2 p d 2 p V(p,p ') exp(ip ■ x — ip ' ■ x '). (A. 7) 

(27r) 4 J 

Introducing 

k = —(p — p'), r = —(x — x"), y = x + x ', (A. 8) 

one obtains 

p-x — p'-x 1 — q-f+k-y, (A. 9) 

such that the momentum-dependent contribution takes the form 

Vp(x, *) = -^- s J^Jdq g - g exp(*g • r) %). (A.10) 

The 5-function arises from the /c-integration and implies x ' = —x as well as f = x, 
which just means that the potential is local in coordinate space. Using 

1 f ,2 gj ~ g| 1 cos(2y) 1 /".a 2gig 2 1 sin(2v?) 

(2vr) 2 J q ql + q 2 2 vr r 2 ' (2vr) 2 J Q q\ + g 2 " tt r 2 ' 1 ] 

with r = r(cos</?,siny2) the g-integration results in 



V p (x,x') 



K 



2 



2np s 



!/ P 2 p2A cos(2y9) sin(2y?) 



= ^cos(2Qp + X )) 

4?rp s r 2 vy; v ; 

In the last step we have introduced P = P(cosx, sinx). 
Similarly, the rest-frame potential takes the form 

V (x,x') = ^ 1 / d 2 g d 2 fc ^ 2glfel 2 2 ^ 2 ) exp(ig ■ r) exp(ifc • (A.13) 

The fc-integration results in the second derivative of a 5-function which again implies 
x ' = —x as well as r = x. Hence, also the rest-frame potential is local and one can 
write 

V (x, x') = V^dy^Siy), (A.14) 
with Vijif) implicitly defined through eq. flA. 13|) . In order to figure out how V (x, x ') 
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acts on a wave function we calculate 

($|y |tf) = J d 2 x d 2 x' (®\x)V (x,x '){x'\V) 

= J d 2 x d 2 x' (<$>\x)V l3 {r)d m d yj 5{y){x '|*> 

= | rf 2 r A ($|| + mA^dyAAv)^ ~ W 

= \j d 2 r V t3 {r)d n d r] m?)(-m) 

= \fd 2 r ($|r) [d ri d r M^)} M*>- ( A -!5) 
It is now straightforward to convince oneself that 

W, ( q = ^-^ + ^ = ^M. (A.16) 
4 J 7rp s r 8 irp s r 4 

Altogether, in coordinate space the resulting potential is hence given by 



K 2 

(r'jf'.\V\r + rJ) = 



2np s 



cos(4v?) P 2 cos (2(<p + x)) 



r 4 2 r 2 



<J(f+-fl) 5(f_-f|). (A. 17) 
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